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ABSTRACT 

We present the mass and X-ray temperature functions derived from a sample of more than 15,000 
galaxy clusters of the MareNostrum Universe cosmological SPH simulations. In these simulations, 
we follow structure formation in a cubic volume of 500/i~^ Mpc on a side assuming cosmological 
parameters consistent with either the first or third year WMAP data and gaussian initial conditions. 
We compare our numerical predictions with the most recent observational estimates of the cluster 
X-ray temperature functions and find that the low normalization cosmological model inferred from 
the 3 year WMAP data results is barely compatible with the present epoch X-ray cluster abundances. 
We can only reconcile the simulations with the observational data if we assume a normalization of the 
Mass- Temperature relation which is a factor of 2.5 — 3 smaller than our non-radiative simulations 
predict. This deviation seems to be too large to be accounted by the effects of star formation or 
cooling in the ICM, not taken into account in these simulations. 
Subject headings: cosmology: theory ~ clusters: general - methods:numerical 



> 

o 

m 

rn 
O 

o 



X 



1. introduction 

Clusters of galaxies are strong X-ray emitters that can 
be observed at large distances using the XMM-Newton 
and Chandra X-ray telescopes. They are excellent cos- 
mological probes that can be used to put strong con- 
straints on the matter density of the universe, (^Im), 
the normalization of primordial density fluctuations (cts) 
and the associated spectral index [n). The number 
of massive clusters in cold dark matter-dominated cos- 
mol ogies is known to b e exponentially dependent on 
(Tg (jSeth fc Tormen"2002') . as has been extensively con- 
firmed by simulations. Therefore, the determination of 
the abundance of massive clusters gives one of the best 
constraints on the normalization of the initial power spec- 
trum of density fluctuations, provided we adopt gaussian 
initial conditions. An independent measurement of the 
cosmological parameters comes from the study of CMB 
anisotro pics. The most rece nt data from 3 year WMAP 
satellite ^Spergel et al.ll2007t l (WMAP3) gives a value for 
(Ts ~ 0.76 ± 0.05 to within Icr error. This is smaller than 
the previous value of = 0.84 ±0.04 estimated from the 
first year WMAP data tSpergel et a l. 2003) (WMAPl). 
This difference in the normalization and the matter con- 
tent — 0.24 vs. VIm = 0.3), translates into large 
differences, up to an order of magnitude as we will show 
in this letter, in the number density of the most mas- 
sive objects formed at present in the Universe. Other 
recent papers independently also ar gue against the low 
values of as obtain ed from WMAP 3 (jEvrard et. al.ll2007t 
iRozo et. al]|2007[) . 



Electronic address: 
Electronic address: 
Electronic address: 



gustavo .yepes@ua iiL.es, raul.sevilla@uam.es| 



sgottloeberSSaip.de 
silkSS astro. ox. ac.ulS 



In order to compare the theoretical cluster mass func- 
tion for a particular cosmological model with the ob- 
served abundance of X-ray clusters as a function of the 
ICM gas temperature, one has to assume that the Mass- 
Temperature relation is sufRciently well known. The 
main obstacle is the accuracy in the determination of this 
relation. Small differences can lead to la rge changes i n 
the determination of cluster mass (see e.g lHenrvl ()2004D ) 
from the X-ray temperature. Estimations of the M-T 
relation from gas dynamical simulations show large dis- 
crepancies, mainly due to numerical r esolution effects 
as we ll as to the physics involved (see lAscasibar et al.l 
(|2006l ) for a review). Most previous numerical studies on 
the comparison of cluster mass functions and X-ray tem- 
perature and/or luminosity functions have either high 
numerical resolution and a low number of objects or 
larger statistics but with very low resolution. 

The aim of this letter is to study the X-ray cluster Tem- 
perature Function (XTF) obtained from a set of large- 
scale non-radiative gas dynamical simulations with suf- 
ficient numerical resolution and statistics to cover the 
range of temperatures for which observational estimates 
of the cluster abundance are known. Our main goal is 
to test whether the observed number of X-ray emitting 
galaxy clusters can be obtained in a cosmological model 
with parameters consistent with WMAP3 or WMAPl 
data at the present time. For this purpose, we compute 
the XTF directly from simulations and compare them 
with the most recent observational estimates. At the 
same time, we derive the values for the normalization of 
the M-T relation that best fit the simulation mass func- 
tions to the observed XTF and compare them with the 
M-T resulting from the simulations. 
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TABLE 1 

Main features of the simulations used in this work, a and logAfo are the best fit parameters of the log(A/2oo/-'*'^o) = 
a log(Tx/3keV) RELATION. Errors correspond to la of the linear fit. AlogAfo is the maximum scatter in the normalization 

(see text) 



Name 


Np 


Hm 


h 


eg 


n 


a 


log Mo 


A log Mo 


MUG 


2 X 1024^ 


0.3 


0.7 


0.9 


1 


1.89 ±0.02 


14.64 ±0.01 


0.46 


MUCL 


2 X 5123 


0.3 


0.7 


0.9 


1 


1.71 ±0.02 


14.52 ±0.01 


0.35 


MUWHS 


2 X 5123 


0.24 


0.73 


0.8 


0.95 


1.62 ±0.05 


14.56 ±0.01 


0.31 


MUW 


2 X 512^ 


0.24 


0.73 


0.75 


0.95 


1.65 ±0.04 


14.54 ±0.01 


0.28 


MU2W 


2 X 512^ 


0.24 


0.73 


0.75 


0.95 


1.65 ±0.05 


14.56 ±0.01 


0.28 


MUW+MU2W 


2 X 512^ 


0.24 


0.73 


0.75 


0.95 


1.60 ±0.04 


14.54 ±0.01 


0.28 



2. SIMULATIONS 

To study the X-ray cluster abundance, we have per- 
formed a series of njMi-radiativc SPH simulations with 
the GADGET2 code (jSpringe l 2005) at the Barcelona Su- 
percomputer Center. Starting at redshift z — 40, we 
followed the non-linear evolution of structures in gas 
and dark matter (DM) to the present epoch {z = 0) 
within a comoving cube of 500ft.~"'^Mpc on a side. The 
so-called MareNostrum universe is the SPH simu- 
lation with 2 X 1024'^ particles (muc ). We assumed a 
concordance cosmological model with the following pa- 
rameters: total matter density flm — 0.3, the baryon 
density fif, = 0.045, cosmological constant Ha = 0.7, 
Hubble parameter h = 0.7, slope of the power spectrum 
n = 1 and normalization erg = 0.9. We also ran the 
same simulation with exactly the same initial data but 
lower mass resoluti o n (2 x 512^, mucl ) as described in 
iGottlober Yepesl (l2?)07h . 

After the release of the 3-year WMAP data, we comple- 
mented our numerical data set with new simulations of 
the same computational box but using WMAP3 cosmo- 
logical parameters: = 0.24, ilt = 0.0418, = 0.76, 
h = 0.73, n = 0.95 and erg = 0.75. We have changed 
both r^jvf and erg (rather than only erg) so as to remain 
on the WMAP degeneracy line for these two parameters. 
As in the concordance model, the power spectrum was 
kindly provided by Wayne Hu, who computed it by direct 
numerical integration of the Boltzmann code. We gen- 
erated the initial conditions for the WMAP3-compatible 
simulations with 2 x 512^ (muw ) particles in exactly 
the same way as for the MareNostrum Universe. In 
order to study the effects of cosmic variance, we have 
completed a second simulation with a different random 
realization (mu2w ). Furthermore, and driven by the 
results obtained for the XTF from these simulations, we 
have also repeated the MareNostrum Universe re- 
alization of the WMAP3 cosmology, but with a higher 
normalization of the initial power spectrum (cig — 0.8), 
consistent within la with the WMAP3 best fit (muwhs ). 
In Table [U we summarize the main characteristics of the 
simulations and the corresponding acronyms for refer- 
ence in what follows. The best fit values of the Mass- 
Temperature relations from clusters obtained in each 
simulation are also shown in the last two columns (see § 
|4]). The clusters have been identified in the simulations 
by means of a hierarchi cal Friends-of-Friend s (EOF ) halo 
finder as described in lGottl5ber fc Yepei ()2007f ). For 
comparison with observational data, we have estimated 
total masses (dark + gas) of clusters at different spher- 
ical overdensities (200, 500, 2500) with respect to the 
critical density. To this end, we started at the position 
of the most massive substructure of the clusters iden- 
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Fig. 1. — Cumulative mass function from the simulations de- 
scribed in the text and Table [T] 

tified with FOF and used th e Bound Density Maxima 
algorithm (|Klvpin et al.lll999| ) to find the spherical over- 
densities. 

3. CLUSTER mass FUNCTIONS 

In Fig. [T] we plot the resulting cumulative mass func- 
tions for all the simulations described in Table [TJ In this 
figure, the total mass of objects corresponds to the re- 
gion enclosing an overdensity of 200 around the center 
of mass found as described in the previous section. As 
can be deduced from this figure, there are no significant 
resolution effects on the number of objects as a function 
of mass. The mass functions for simulations MUC and 
MUCL nicely overlap each other despite the fact that they 
differ by a factor of 8 in mass resolution and a factor of 
^ 3 in spatial resolution. On the other hand, there is 
a significant difference in the number of cluster-size ob- 
jects depending on the cosmological model. The number 
density of clusters with masses M200 > 5 x lO^^/i~^M0 
in both simulations with the low normalization, best-fit 
WMAP3 cosmological parameters, MUW and Mu2w , is 
^10 times smaller than for the simulations of the con- 
cordance ACDM model. The muwhs simulation with 
fjg = 0.8 has a number density 2 higher than the 
simulations with tig = 0.75, but still is a factor of ^ 5 
smaller than in the concordance cosmology. Finally, Fig 
[1] also shows that the effects of cosmic variance is not 
important in determining the abundance of clusters at 
these scales. The agreement of the mass functions for 
the two different realizations of the WMAP3 cosniologi- 
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cal model clearly confirms this. On the other hand, we 
have also checked for the possible effects of small volume 
sampling in the determination of the mass function for 
the most massive objects. To this end, we have com- 
pared our mass function for the muwhs simulation with 
the mass function obtained from a dark matter only sim- 
ulation of the same cosmological model as MUWHS and 
number of particles but larger computational volume (1.5 
Gpc). This simulation has been done also at MareNos- 
trum with the GADGET2 code by P. Fosalba for the 
Dark Energy Survey project. The agreement between the 
two mass functions is remarkable for halos with masses 
M200 > 5 X IO^^H-^Mq. 

Therefore, we conclude that the estimation of the clus- 
ter mass function from our simulations is robust and not 
likely to be affected by numerical effects. Now, our pur- 
pose is to compare them with data coming from X-ray 
observations of clusters. As our simulations include gas 
dynamics, we can directly measure the X-ray temper- 
ature from the gas content of our halos. In DM-only 
simulations, one has to rely on the Mass- X-ray Temper- 
ature relation to transform mass into temperature or vice 
versa. Here we will do the same exercise and compare 
the calculated XTFs. 

4. X-RAY TEMPERATURE FUNCTION 

The most recent published data for the XTF of nearby 
clusters uses temperatures derived fr om X-ray observa- 
tions i nainly by the ASCA satellite (jlkebe et al.l 120021: 
lHenrvl[200l as a measure of the mean temperature of 
the ICM. The differences shown in the temperatures of 
clusters from these two datasets reflect the systematic er- 
rors in the observed XTF. For our simulated clusters, we 
computed several temperature estimations: the emission- 
weighted temperature, T^w by weighting the tempera- 
ture of each SPH particle within the cluster by their 
X-Ray luminosity. We also computed the spectroscopic 
temperature, Ty , following the procedure described in 
IVikhlininI (|2006l ). which is suppose to give a more ac- 
curate value of the observed temperature of an X-ray 
emmitting plasma. Therefore, in what follows we will 
use Tx for the simulated clusters. 

In Fig. [2] we show the cumulative XTF as a function 
of the spectroscopic Tx for the clusters found in simu- 
lations described in Table [H We also represent the ob- 
servational d ata as points with error bars as described in 
iHenrvl (|2004D . The observational data were rescaled to 
units of /i = 1. 

The predicted number density of X-ray clusters above 
a given temperature for the MUC and mucl simulations 
with cTg — 0.9 is in good agreement with the data. Again, 
as in the case of mass, the WMAP3 most favored cosmo- 
logical model underpredicts the density of X-ray clusters 
with respect to the observations by a factor of ~ 10 for 
clusters with > A keV. The situation is slightly bet- 
ter for the higher normalization muwhs simulation. But 
still, it predicts a factor of ~ 6 fewer density of clusters 
hotter than T^^ > 4 keV than in reality. 

We showed in Fig [1] that effects of resolution are neg- 
ligible in the estimate of the cumulative mass function 
for massive clusters. This could not be the case for the 
temperature estimates from the gas particles. In order 
to check whether the XTF could be affected by resolu- 
tion, we also show in Fig [2] a comparison of the XTF 
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Fig. 2. — Cumulative X-Ray Temperature function for two dif- 
ferent cosmologies. Same notation as in Fie. Ill Solid triangles 
and squares represent Ikcbc ct al. (2003) and IHenrv (2004) data, 
respectively. Errors were computed as explained in .Henrvi t2004i) . 
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Fig. 3.— Best-fit cumulative XTF to data from WM APS simula- 
tion mass functions (see Fig[T]l and varying the two free parameters 
of the Af200 ~ relation. The XTF from the MUG simulation is 
shown for comparison. 

between MUC and mucl simulations. As can be seen, 
the spectroscopic temperature estimate of clusters is bi- 
ased high when low mass resolution is used in a SPH 
simulation. Thus, we expect that the difference in XTF 
shown between MUC concordance model simulation with 
1024"^ particles and the WMAP3 lower resolution simu- 
lations (512'^) is in fact a lower limit. If we increased the 
mass resolution of the latter, we would obtain a larger 
difference with respect to the MUC and data. 

5. DISCUSSION 

We have shown in the previous sections that the low 
normalized WMAP3 cosmological simulations underpre- 
dict the abundance of X-ray clusters by a factor that 
ranges between 6 — 10 with respect to estimates from 
ASCA observations. Now, our estimates are based on 
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the results from non-radiative gas dynamical simulations 
of the ICM. There is still no clear answer to what extent 
cooling and star formation are important in the thermo- 
dynamics of the ICM. The extreme complexity of the 
processes involved presents a serious challenge for sim- 
ulating them accurately in a cosmological setting. Re- 
sults from simulations which incorporate some model- 
ing of these processes have shown that the M-T rela- 
tion is not strongly affec t ed by non-gravitational heat- 
ing (|Borgani et aLl l2004t |Nagai, Vikhlinin & Kravtsov 
|2007() . A rather more important ingredient in the deter- 
mination of the XTF from mass functions is the intrinsic 
scatter of the M-T relation. If the scatter is big, then a 
rather low normalization power spectrum can in princi- 
ple give high enough XTF to be compatible with observa- 
tions. Given the very good statistical sample of objects in 
our simulations, we can reliably estimate not only the M- 
T relation but also the intrinsic scatter due to the cluster 
dynamics. In Table [1] we report the least-square fit val- 
ues of the Maoo/Mo/i-i = (Mo/MQ/i-i)(rx/3keV)" for 
the different simulations where errors in both the slope 
a and normalization Mq correspond to la in the fit. We 
can also make a reliable estimate of the intrinsic scatter 
in the M200 — Tx relation. The linear fit of the log Af2oo 
versus logTx has a Pearson's correlation coefficient bet- 
ter than 0.99 for all simulations. The maximum intrinsic 
scatter, AlogMg, is also shown in Table [TJ It is defined 
as the value for which 99% of all the clusters used in the 
fit have their spectroscopic temperature within the values 
logTx = (logM2oo)/a - logMo T AlogMo. As can be 
seen, the values of the scatter are between 0.28 — 0.31 dex 
(~ 2 factor with respect to Mq) for the WMAP3 simula- 
tions. But are the differences shown in Fig[2] between the 
simulated XTF and data compatible with this intrinsic 
scatter of the M-T relation? In order to give a possible 
answer to this question, we have estimated the Mo and a 
parameters of the M200 — Tx relation needed to accom- 
modate the mass functions shown in Fig[l]to the observa- 
tional XTF data by a minimization. We show in Fig[3] 
the best-fit simulated XTF for the WMAP3 simulations 
to the observational data points, together with the simu- 
lation results for the high-normalization MUC simulation. 
The best fit values found for the MUW -I- Mu2w are 
a' = 1.64 and logA/^ = 14.09 for the Ikebe et al data 
and a'^ = 1.49 and logA/g = 14.17 for the Henry data. 
When both observational data sets are taken together in 
the fit, we obtain a^+'' = 1.64, logAfo^+'' = 14.10. For 
the higher normalization WMAP3 simulation muwhs , 
we find = 1.66; logAf^ = 14.18 for Ikebe, a'' = 
1.44 and logA^o'' = 14.28 for Henry and a'+'' = 1.67, 
logM^^^ = 14.18 for the combined datasets. Now, if 
we fix the slope, a, to the best fit value obtained from 
each simulation (see Table [Ij we find a value for normal- 



ization parameter logAfo = 14.10 - 14.13 for WMAP3 
fjg = 0.75 simulations and logA/o = 14.19 - 14.22 for 
the as = 0.8 WMAP3 simulation. Finally, if we as- 
sume the self-similar behavior of the M-T scaling rela- 
tion, a = 3/2, then, the best fit values for A/q are quite 
similar: logA//o = 14.20 - 14.26 for the muwhs sim- 
ulation and logAfo = 14.12 - 14.17 for the MUW + 
Mu2w simulations. Therefore, the normalization of the 
A^200 ^ Tx relation needed to fit the observational XTF 
for the as = 0.75 MUW and Mu2w simulations is a fac- 
tor of 0.40 - 0.45 dex (~ 2.5 - 2.8 times ) smaller than 
the best fit values shown in Table [T] For the as = 0.8 
WMAP3 MUWHS simulation this factor is 0.36 - 0.39 dex 
(r-^ 2.3 — 2.4 factor). As we have seen, the maximum 
scatter derived from our WMAP3 non-radiative gas dy- 
namical simulations is ~ ±0.28 — 0.31 dex (i.e. a factor of 
'--^ 2 ) . It is not clear that non-gravitational heating could 
affect the thermodynamics of the ICM in such a way that 
this could account for a factor of ~ 2.5 — 2.8 lower normal- 
ization with respect to the predictions of the simulations 
reported here. For instance, the normalization for the 
emission- weighted A^soo — Tew frort i SPH simulations in - 
cluding cooling and star formation iBorgani et al.l ()2004l ) 
is a factor of 1.46 smaller than the value we obtained for 
our MUWHS simulation. If we compare the normalization 
of the spectros copic A^son — from the radiative c luster 
simulations of iNagai. Vikhlinin fc Kravtsov! (|2007l ) with 
ours, the difference is within a factor of 1.5 — 1.6. 

In conclusion, it seems unlikely that we can repro- 
duce the observational estimates of abundance of X- 
ray clusters with a normalization of the power spec- 
trum as low as the best fit value given by WMAP3. A 
slightly higher normalization of trg = 0.8 alleviates the 
problem, although the cluster abundance still lies be- 
low the observational estimates. Considerably steeper 
slopes and lower normalization of the M-T relation are 
needed to reconcile the predicted mass functions of clus- 
ters with the observed XTF in this case. Alterna- 
tive explanations which retain a low normalisation of 
as appeal to the effects of primordial non-gaussianity 
(Sadch, Rephaeli & Silk 20 0^ or to dyn amical dark en- 
ergy (Bartelmamij.^oraiij^_^ttinc3|2^ However for 
the standard cosmological model. X-ray clusters of galax- 
ies seem to prefer a higher erg than predicted by the CMB 
anisotropics, in agreem ent with the abu ndance of optical 
clusters from SDSS (Rozo et. al.ll2007D . 
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